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Abstract. The matter content of relativistic jets in AGNs is dominated by a mixture of protons, 
electrons, and positrons. During dissipative events these particles tap a significant portion of the 
internal and/or kinetic energy of the jet and convert it into electromagnetic radiation. While 
leptons - even those with only mildly relativistic energies - can radiate efficiently, protons need to 
be accelerated up to energies exceeding 10 16-19 eV to dissipate radiatively a significant amount 
of energy via either trigerring pair cascades or direct synchrotron emission. Here I review various 
constraints imposed on the role of hadronic non-adaiabatic cooling processes in shaping the high 
energy spectra of blazars. It will be argued that protons, despite being efficiently accelerated 
and presumably playing a crucial role in jet dynamics and dissipation of the jet kinetic energy 
to the internal energy of electrons and positrons, are more likely to remain radiatively passive 
in AGN jets. 
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1. Introduction 

Production of relativistic jets in AGN is very likely mediated by rotation of large scale 
magnetic fields in magnetosphere of black hole and/or accretion disk (Blandford 1976; 
Lovelace 1976; Phinney 1983; Camenzind 1986; McKinney & Blandford 2009). This leads 
to a formation of Poynting flux dominated outflows. Those in turn can at some point 
be converted to the matter dominated jets (Komissarov et al. 2007; Tchekhovskoy et al. 
2009; Lyubarsky 2010; Komissarov 2010), with a terminal Lorentz factor T ~ Pj/Mc 2 
where Pj is the rate of energy extraction from rotating BH and/or accretion disk and M 
is the mass flux. Depending on whether a jet is launched in the BH magnetosphere or by 
the accretion disk, the mass flux is expected to be dominated by electron/positron pairs 
or by protons. 

Presence of protons in AGN jets is indicated by the low-energy cutoffs in the radio 
spectra of hot spots in radio-lobes (Blundell et al. 2006; Stawarz et al. 2007; Godfrey et 
al. 2009) and by circular polarization and Faraday Rotation in radio cores (Vitrishchak 
et al. 2008; Park & Blackmail 2010). If so, such protons might manifest their presence 
also via radiative contribution to high energy spectra in blazars, by synchrotron emission 
of pair cascades triggered by photo-meson process and by direct synchrotron emission of 
protons and mesons (Mannheim & Biermann 1992; Mannheim 1993; Rachen & Meszaros 
1998; Aharonian 2000). Whether such a contribution can be significant is the main is- 
sue of this presentation. We start with a short review of the leading hadronic models 
(§2); discuss efficiencies of the proton acceleration and cooling processes (§3); present 
observational constraints on hadronic contribution to blazar spectra (§4). And, finally, 
we discuss the role of protons they play in the leptonic models regarding the aspect 
of preheating electrons up to thermal proton energies required to participate in further 
stochastic acceleration process (§5). Main results are summarized in §6. 
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2. Hadronic models - basic features 

2.1. Luminous blazars - 'photo-meson' models 

Luminous blazars are hosted by quasars which, when observed away from the axis of the 
jet, form population of FRII type double radio sources. They are powered by accretion 
onto BHs with masses typically of the order of 10 9 solar masses and their accretion 
luminosities are in the range 10 46 — 10 47 erg s _1 . However, when such jets are oriented 
close to the line of sight, these quasars are seen as being dominated by nonthermal 
radiation of jets, with the apparent luminosities 10 48-49 erg s _1 and spectra dominated 
by the broad high energy component peaking in the 1-100 MeV band. This radiation 
shows high amplitude variability on time scales from years down to days and even hours, 
implying the strong dissipative events taking place not far from the base of a jet. There the 
energy densities of the magnetic and radiation fields are very large, providing conditions 
for acceleration of protons up to energies of 10 9-10 GeV, and for cooling them via inelastic 
collisions with soft photons (Sikora et al. 1987). For typical background radiation fields 
most of proton energy is converted to mesons and this initiates processes which according 
to proposers of hadronic models are responsible for 7-ray production in luminous blazars 
(Mannheim & Biermann 1992). 

These processes are dominated by the following channels. In approximately 90% col- 
lisions, the produced mesons are pions. In 2/3 of them they are the neural pions (ir°) 
and in 1/3 of them - the positive pions ir + . Neutral pions almost immediately decay into 
photons which in turn trigger pair-cascades driven by photon-photon pair production and 
their synchrotron radiation. Escaping radiation is the product of 3rd and 4th generation 
of pairs (radiation of the first two are totally converted to e+e - -pairs). The resulting 
electromagnetic spectrum is predicted to form the high energy component peaked in the 
7-ray band, with a high energy break at hvb r ^10 — 30 GeV, where r 77 ^ e + e - ~ 1, and 
low energy tail in the X-ray band with a slope ax > 0.5. 

Such spectra can be affected, but not significantly, by electromagnetic output resulting 
from production of charged pions via processes pj — > n + tt + and nj — > p + 7r~. The 
charged pions decay producing muons (7r ± — > /i ± + v^jvp), and the muons decay pro- 
ducing positrons/electrons — > e^ + D e /v e + v il jv ll ). The resulting electrons/positrons 
join the pair cascade triggered by the decays of the neutral pions and together with 
the synchrotron emission of muons (Rachen & Meszaros 1998) contribute an additional 
~ 6-16 % to the electromagnetic output of the photo-meson process. 

2.2. Low luminosity TeV BL Lac objects - proton- synchrotron models 

The low luminosity BL Lac objects arc hosted by radio-galaxies of type FR I. Just as 
FR II radio sources, they are associated with giant elliptical galaxies, with central BH 
masses of the order of 10 9 M Q , but with the jet powers at least 3 orders of magnitude 
lower than in the radio-loud quasars and extremely low accretion luminosities. Doppler 
boosted nonthermal radiation from their jets reaches TeV energies, and the luminosity 
peak of the high energy spectral component is located at GeV energies. Resulting from 
the low radiative environment, energy losses of ultra-relativistic protons in theses objects 
are presumably dominated by direct synchrotron emission of protons. This mechanism 
was suggested by Aharonian (2000) to be the primary source of 7-rays in low luminosity 
BL Lac objects. Main spectral features of such models, as is in the case of synchrotron ra- 
diation of electrons, are directly related to the magnetic field intensities and the injection 
function of rclativistic particles (here of protons). They produce photons with average 
energies v p ^ syn — (2e/ (3irm p c)) j p B"D and energy spectra with the slopes a — (q p — l)/2 
in the slow cooling regime and a = q p /2 in the fast cooling regime, where B' is the 
intensity of magnetic field in the jet co-moving frame, V is the Doppler factor, q p is the 
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index of the proton power-law injection function, Q p oc j p " , and a is the index of the 
radiation energy flux, F v oc v~ a . Application of those formulae indicates that protons 
accelerated up to Lorentz factors j p ~ 10 10 and cooled in magnetic fields B' <~ 100 Gauss 
may produce spectra reaching TeV energies and with the luminosity peak located near 
Vp,s V n,max- An investigation whether these parameters are feasible - given the constraints 
imposed on -f p , m ax by jets with the limited power and magnetisation - is given in §4.2. 



3. Radiative efficiencies 

3.1. Time scales 

Acceleration . 

Time scale of the proton acceleration as measured in the co-moving frame of the flow 

is 

(cc = We=^^, (3-1) 

where Rl is the Larmor radius, B' is the magnetic field intensity, and / is the parameter 
which in the case of shock acceleration depends on the spectrum of magnetic turbulence 
and on the velocity of the upstream-flow (Rieger et al. 2007) and for mildly rclativistic 
shocks is expected to be at least of the order of 10. 
Adiabatic losses 

Assuming that cross-sectional radius of the source is of the order of the cross-sectional 
radius of the jet, R, relativistic plasma moving down the jet with a Lorentz factor T 
undergoes energy losses due to adiabatic expansion. For conical jets with the half-opening 
angle 0j = R/r < 1/T, where r is the distance of the source in a jet, the time scale of 
the adiabatic losses is (Moderski ct al. 2003) 

, 3 R . 

tad = 2 Te~ry c ■ (3 > 

Photo-meson process 

Time scale of the energy losses via the photo-meson process can be estimated using 
the approximate formula (Begelman et al. 1990) 

t>P1 ~ (° P -,K p ,)cn' ph { V > > v' th) ' (3 - 3) 

where {a pl K pl ) ~ 0.7 x 10~ 28 cm 2 is the product of the photo-meson cross-section 
and inelasticity parameter averaged over the resonant energy range, n' ph {v' > v' th ) = 
J , n! v dn' v , hv' th ~ m 7r c 2 /7 p is the threshold photon energy and is the rest mass of 

th 

the pion. The target radiation field is provided by the internal and external sources. The 
internal seed soft radiation is dominated by synchrotron emission of primary (directly 
accelerated) electrons, the external one - by re-scattered/reprocessed disk radiation. Ap- 
proximating the broad synchrotron spectral component by a power-law function with the 
energy- flux index a = 1 and denoting its luminosity by L s , we have 

Spectra of external radiation fields are in turn narrow and therefore can be approxi- 
mated by mono-energetic functions. Hence, 

£,L d T JL^Tf 2 

~ W h Vext - A,cR^h Vext f ° r > > m « C ^ ' (3 - 5) 
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and n' p h t ext = f° r 7p < m^c 2 /(Thv ex t), where is the luminosity of the accretion 
disk and £ is the fraction of this luminosity rescattered/reprocesssed on a spatial scale 
corresponding with a distance r of the source in a jet. 
Synchrotron emission 

Time scale of the proton cooling via the synchrotron process is 

, 3 /t7j p \ 3 m e c 1 

where u' B = B' 2 /(8ir) is the magnetic energy density. 

3.2. 'Mono-energetic' efficiencies 
In order to illustrate efficiencies of the proton acceleration and cooling processes we 
introduce their dimensionless rates, as scaled by the adiabatic losses rates, tj = t' i 1 /t' a d- 
They are: 

_ 7 B'R {a/Q.lf/ 2 L)% l 

Tacc " 48x10 mrW P " °' 8 (//io) (r/io) fcf) ^ ■ (3J) 



L2 x 10 (^r) v*R~ u (OjT) (r/io) 4 i? 16 7p4 ° ' {6± 

Jext) 2Sx1n -40 iMgjT) ^ (^)4 5 (g 3 r) 

P7 ~ ' ^ t ri? " ' (/uWioeVHr/io)^' 



(3.9) 



and 



_ 29 RB' 2 lp nQ (<7/0.1)L jl4 7 



' p,syn 



1.1 x 10-- c 0.3 - ^ ; 7p , 10 , (3.10) 



r " " (^r)(r/io) 2 i?i 6 



where Lb — cu' b ttR 2 T 2 = aLj/(l + a) is the magnetic energy flux, Lj = Lm + Lb is 
the total jet power, a = Lb/Lm, and Lm is the energy flux of the rest mass. Assuming 
a < 1 we use approximation Lb ~ er£ j • 

The 'scaled' quantities, £7,47, 7 P .io, -R16, -^s,47, and {£,Ld)ia are defined in the usual 
way, i.e. X n = X/10 n . 

3.3. Maximal proton energies 
Maximal proton energies - if limited only by adiabatic losses - can be found from 
T acc = 1 to be 

q (a/O.l) 1 /^ 1 / 2 
7,(^ = 1)^8x10° (//1Q)(r/10) 3 ( g r) ' 

Stronger limits are imposed if dominant energy losses are non-adiabatic, i.e. for t coo ; = 
t P7 + Tp^ S y n > 1. For energy losses dominated by photo-meson process with the target 
radiation field provided by internal sources or by proton-synchrotron radiation 



7 P ,max = l P {T acc = 1) Min[l; ^R/Rc} . (3.12) 
where in the 1st case (photo-meson process) 



and in the 2nd case (proton-synchrotron radiation) 



(3.13) 



1« W0-l) 3/2 ^47 

^^"^ (/W^fffl^ - (3 ' 14) 
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For energy losses dominated by photo-meson process with the target radiation field pro- 
vided by external sources, T p ^ xt>> is predicted to be lower than unity for any proton energy 
and therefore jp^max is not expected to be affected. 

3.4. Total efficiencies 
Radiative efficiency of a given cooling process can be estimated using formula 

f^mhp^ipdip 

Ji Q~i P lpdl P 
where rji(j p ) ~ Min[Ti; 1] and Q 7p is the proton injection function. 

We calculate such efficiencies below assuming power-law injection of protons Q 7p oc 
7 P 9p with q p — 2. This specific value of the index is chosen because efficiencies obtained 
for q p — 2 provide upper limits of efficiencies available for q p > 2 being predicted by 
theoretical models of a diffusive acceleration of particles in mildly relativistic shocks and 
supported by the slopes of synchrotron spectra produced in the IR-Optical bands by 
primary (directly) accelerated electrons/positrons. 

For energy losses dominated by photo-meson process with internally produced seed 
photons and by proton-synchrotron emission, we obtain 

Vi ~ 1 + for R < fl c; (3.i 6 ) 

m "fp,max 

and r\i < 1 / \n-y p ^ max for R> R Cl where in case of 'internal' photo- meson process, R c is 
given by Eq.(13) and 

7 ™ = 7p ( T #*> = 1) * 8.3 x 10 s (^ r ))(r/10) 4 fi 16 j (3 1?) 

while in case of proton-synchrotron emission, R c is given by Eq.(14) and 

H Vn) = 7,(w) - 3.4 x 10- Vf^f* 1 * ■ (3.18) 

For energy losses dominated by photo-meson process with externally produced seed 
photons, at any distance larger than 



YR^f = 1) 1 inl6 (hv ext /WeV) (r/10) 2 

(ext) 



( P>7 < 1 and the efficiency is 

^ext) ^ ^(ext) m {ip.max/lp.th 

where lp , th = m«c 2 / {Thv ext ) ~ 1.4 x 10 6 /((r/10)(/«Wl0cV)). 



„{ext) ^ Aext) 111 V Ip.max/ Ip.th) ,„ „n 

Vpi - T p~i l„ ' (J.ZUj 



4. Observational constraints 

4.1. Luminous blazars 

In order to account for the 7-ray luminosities of powerful blazars, radiative efficiency 
has to be 

^ >0 - 3 (r/io)MWo.3)W (41) 
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where r/diss is the fraction of the jet energy flux dissipated in the 'blazar zone'. As it can 
be verified using approximate formulae presented in §3.4, such efficiency is difficult to 
reach even for protons injected with the energy distribution slopes q p — 2. For our fiducial 
parameters it is about 10 % for the photo-meson process with intenally produced seed 
photons and much less for others. For slope q p ~ 2.4 the efficiency drops to r\ ~ 10~ 3 . 

Hadronic models may have also problems to explain very hard X-ray spectra of lu- 
minous blazars. Those blazars often have slopes a x < 0.5 (see Table 1 in Sikora et al. 
2009) and in order to explain such spectra by synchrotron radiation of secondary - 
products of the cascades powered by hadrons - one needs to assume inefficient cooling 
of ultra-relativistic electrons/positrons up to energies 

n « ^(WfOOkeV) 1 / 2 , ln . 

Inefficient cooling of such energetic electrons/positrons implies very weak magnetic holds 
and, therefore, puts strong constraints on the efficiency of the proton acceleration and 
on efficiency of photo-meson energy losses via limitation of the maximal proton energy. 
Furthermore, as Sikora et al. (2009) demonstrated, in order to avoid overproduction of 
X-rays in these magnetically weak sources by SSC radiation of primary electrons, it is 
necessary to assume the source sizes of the order of parsecs, and for such sources the 
efficiency of the photo-meson process is further reduced. 

4.2. Low luminosity BL Lac objects 

In these objects, because of low radiation energy densities - both in jets themselves 
and in the surroundings of the jet - the non-adiabatic energy losses of protons are pre- 
sumably dominated by the proton-synchrotron mechanism. However, noting that such 
objects are hosted by weak radio-galaxies, with the jet powers Lj ^ 10 44 crg s _1 , ef- 
ficiency of synchrotron-proton models is also expected to be strongly reduced because 
weaker magnetic fields. In order to keep them at reasonable level more compact sources 
must be assumed. However, even in the case of most relativistic protons, the required 
size of the source to provide sufficiently strong magnetic fields for efficient cooling is 
unreasonably small (see Eq. (3.14), 

R ^°- 7><io " wMwnli)3 '"''- <«> 

This is 3 orders less than the gravitational radius of the BH with mass Mbh ~ 1O 9 M0. 
Considering the minimal cross-sectional radius of a jet to be R ~ 10 15 cm, which for 
8j ~ 1/r corresponds with a distance 100 gravitational radii of the 10 9 M Q BH - required 
to be at least of this order to accelerate the jet up to T ~ 10 (Komissarov et al. 2007) - 
we can find using Eqs. (3.10) and (3.11) that 



(q/0.1)£ 3 -4 4 ) 3 

(//i)(^r) (r/io^ 



T p ,syn{lp, m a X ) 0.7 X 10" 4 - ; ,\ '- Vwwi ^3 o ' ( 44 ) 



This indicates that even for such extreme parameters as / ~ 1 and a ~ 1 the efficiency is 
too small to explain 7-ray luminosities L 1 ~ 10 erg s _1 , unless one assumes very hard 
(q p < 1) proton injection spectra and adopts significantly larger total jet power. 

The main purpose of proton-synchrotron models was to explain the relatively stable 
shape of the TeV spectra in variable low luminosity BL Lac objects (Aharonian 2000). 
Obviously, the critical issue of such models is whether protons can reach sufficiently 
large energies to produce synchrotron spectra extending up to TeV energies. Combining 
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formulae for average synchrotron photon energies, magnetic energy flux, and maximal 
proton energies (Eq. 3.11) gives 



(a/0A)^L%\ 
(//10) 2 (F/10) 2 (tf,F) 2 i? 15 



hlSp, S yn,max - 3 X 10 5 fr/lCW? fflW d I TcV ]- ( 4 - 5 ) 



One can see that spectra may extend to TeV energies only if assuming / ~ 1, a ~ 1, and 
jet powers Lj > 10 45 erg s^ 1 . 



5. Protons in leptonic models 

Disproving interpretation of high energy spectra of blazars produced via hadronic 
models does not disprove the presence of protons in AGN jets. They simply are expected 
to be radiatively inefficient but are likely to dominate the jet energy flux and strongly 
affect the dynamics of dissipation processes via shocks in the regions of low magnetization 
parameter (er < 0.1). In this scenario, a large fraction of dissipated energy must be 
converted to relativistic electrons, otherwise leptonic models will be inefficient despite 
high radiative efficiencies of energetic electrons. However in order to get electrons to 
participate together with protons in the diffusive shock acceleration process, electrons 
must be first preheated up to average energy (to be strict - up to average momentum) 
of protons heated by randomization and compression in the shocked plasma. Several 
scenarios have been investigated for such electron preheating, both analytically and in 
PIC-simulations (Amato & Arons 2006; Amano & Hoshino 2009; Sironi & Spitkovsky 
2010). Some of them indicate the formation of a power-law energy distribution with 
the slope 1 < q e < 2. If it is true, then the number of electrons joining protons in the 
stochastic acceleration process can be lower than the number of protons by significant 
factor. Hence, the high efficiency of blazar radiation may indicate that q e < 1, or that 
there is a significant pair content. 

In the ERC (Extcrnal-Radiation-Compton) models (Sikora et al. 1994) of 7-ray pro- 
duction in luminous blazars such a low-energy tail should be observed in the 30keV-l 
MeV spectral energy range. Unfortunately, these bands are observationally poorly cov- 
ered, particularly above 20 keV. A number of blazars have been detected up to ~ 50keV, 
by INTEGRAL (Beckmann et al. 2009), by Swift/BAT (Ghisellini et al. 2010), and by 
others (see Table 1 in Sikora et al. 2007). Some of them have spectral slopes correspond- 
ing with 1 < q e < 2. This, together with the location of the break in the energy range 
1 - lOMeV as suggested by CGRO/COMPTEL and OSSE observations (Zhang et al. 
2005; McNaron-Brown et al. 1995), may indicate a moderate (n e /n p ~ 10) pair content, 
implied also by studies of the bulk-Compton and Compton-rocket effects (Sikora & Made- 
jski 2000; Ghisellini & Tavecchio 2010). At energies < 20 keV most blazars have softer 
X-ray spectra, but they can result from the contribution of the SSC process and/or from 
the superposition of X-rays produced at different locations in a jet, such as the orphan 
X-ray outburst detected in 3C279 (Abdo et al. 2010). 



6. Conclusions 

For realistic jet powers (Lj ^ LEdd in luminous blazars and Lj ^ 10~ 3 Lsdd in low 
luminosity BL Lac objects), limited magnetization (a ^ 0.1 is required to allow formation 
of strong shocks), and moderately steep proton injection function (q p ^ 2 is suggested 
by IR-optical spectra and by theoretical acceleration models), the hadronic models fail 
to: 
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• reproduce 7-ray luminosities of blazars; 

• explain formation of very hard X-ray spectra in luminous blazars; 

• provide the spectral extension up to TeV energies in low luminosity blazars. 

Nevertheless, as indicated by several independent observations, protons are present in 
AGN jets and presumably play key role in dissipation processes in shocks. In particular, 
they transfer a fraction of the dissipated energy to electrons/positrons helping them to 
reach threshold energies for further acceleration by stochastic mechanisms and produce 
the observed 7-ray spectra via the ERC and SSC scenarios. However, significant pair 
content may be required to achieve reasonable efhecincy of that energy transfer. 
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Discussion 

BEDNAREK: Can curvature energy losses of protons be important in AGNs? 
SIKORA: They cannot be. This is because following particle acceleration processes in 
shocks or reconncction regions, relativistic protons are injected with a broad distribu- 
tion of pitch angles. Such protons (as well as all other charged particles) spiral around 
magnetic field lines rather than slide on them. 

YUAN: There are two kinds of B field in a jet. One is a large scale helical field, another 
being turbulent in a shock front. When we calculate radiation spectrum, we use turbulent 
one, but when we calculate polarization, we seem to use the large scale ordered field. 

SIKORA: Properly calculated polarization must take into account both, the tangled/turbulent 
magnetic fields compressed and/or generated in a shock, as well, as the large scale mag- 
netic fields transmitted through the shock. 

DERMER: Acceleration to high energies is faster for ions than protons. How do your 
conclusions change if you consider Fe rather than p? 

SIKORA: For approximately solar abundances of AGN plasmas the number of heavy 
nuclei is 3> Z times smaller than of protons. Therefore, despite the fact that heavy 
nuclei are accelerated Z times faster, they contribute to radiative processes much less 
than protons and our main conclusion that hadronic models cannot reproduce 7-ray 
luminosities of blazars remains valid. 

PIRAN: Auger indicates that UHECRs are nuclei. It is much easier to calculate nuclei. 
Nuclei will diffuse in the intergalactic field and can propagate only up to ~ 10 Mpc. 
This suggestes that Cen A is the main source of UHECRs if those are nuclei and if the 
intergalactic magnetic field is > 1CT 9 G (Piran et al. 2010). 

SIKORA: Yes, I agree that Cen A may contribute significantly to the observed UHECRs. 



